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We numerically investigate quantum quenches of a nonintegrable hard-core Bose-Hubbard model 
to test the accuracy of the microcanonical ensemble in small isolated quantum systems. We show 
that, in a certain range of system size, the accuracy increases with the dimension of the Hilbert 
space D as 1/D. We ascribe this rapid improvement to the absence of correlations between many- 
body energy eigenstates. Outside of that range, the accuracy is found to scale as either I/'/D or 
algebraically with the system size. 
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Introduction .— The microcanonical ensemble (MCE) 
is the fundamental working hypothesis in statistical me¬ 
chanics for the description of equilibrium states in an 
isolated system [l[. However, the range of applicability 
and the accuracy of the MCE have yet to be fully un¬ 
derstood. This OToblem dates back to von Neumann’s 
seminal work and has seen resurgence of interest 

recently partly because isolated quantum systems 

have been realized using ultracold atoms and thus 

the verification of the hypothesis has now become an is¬ 
sue of practical importance as well. 


The understanding of why the MCE can describe equi¬ 
librium states in the thermodynamic limit (TDL) has 
deepened considerably in recent years. Even under uni¬ 
tary time evolution, an effective stationary state can ap¬ 
pear due to dephasing between many-body energy eigen¬ 
states (Ill - [l3| : then, physical quantities are obtained by 
the weighted average of expectation values over individ¬ 
ual energy eigenstates (see Eq. ([3])) with the weights de¬ 
termined by the initial condition. On the other hand, it 
is empirically known that the MCE, which presupposes 
the equal-weighted average of those expectation values 
(see Eq. (H))), well describes the physical quantities in 
the effective stationary state, regardless of the initial con¬ 
dition. Recent studies have suggested that this success 
of the MCE derives from the fact that quantum states 
with close eigenenergies typically emulate the same ther¬ 
mal state in the TDL. This scenario has been formulated 
mathematically [2, [l5l| and verified by uniform random 

samplings of states in a narrow energy shell fl6l-[l^. The 
eigenstate thermalization hypothesis (ETH) is a 

quintessential representation of this scenario. The ETH 
states that expectation values of a macroscopic observ¬ 
able for energy eigenstates in a narrow energy shell, which 
we shall refer to as eigenstate expectation values (EEVs), 
are equal in the TDL. The ETH has been supported by 
several numerical studies in various nonintegrable mod¬ 
els 


While these studies concern why the MCE works in the 
TDL, we address the question of how accurately it works 


in small isolated quantum systems. Such a question can 
now be tested experimentally [Sll - I^ . We find that the 
ETH gives, in small isolated quantum systems, an up¬ 
per bound on the accuracy of the MCE and the bound 
scales as 1 / y/P, where D is the dimension of the Hilbert 
space [^. However, the accuracy itself has not yet 
been studied in concrete physical systems. 

In this paper, we numerically investigate the accuracy 
of the MCE for quantum quenches in a one-dimensional 
nonintegrable tight-binding model of hard-core bosons 
(HCBs) with the number of sites L = 15,18, 21 and 24. 
We show that, for some quench magnitudes, the accu¬ 
racy improves proportionally to 1/D, which is much bet¬ 
ter than the upper bound given by the ETH mentioned 
above. We argue that the 1/D scaling implies that for 
quantum quenches in nonintegrable systems, no correla¬ 
tions arise between many-body eigenstates. We also show 
that as we increase the system size with the quench mag¬ 
nitude fixed, three distinct regimes emerge in which the 
accuracy scales as (i) I/'/D (ETH regime), (ii) 1/D (no¬ 
correlation regime), and (iii) with a > 0 (algebraic 
regime). 

Formulation of the problem. — We begin by formulat¬ 
ing the problem in a general setup. We consider an iso¬ 
lated quantum system described by a time-independent 
Hamiltonian H. Let {\En)}n=i be the set of the eigen¬ 
states of H with eigenenergies iL|D„) = 

En \En) for n = 1, 2,..., D, where D(^ 1) is the dimen¬ 
sion of the Hilbert space. An initial pure state I'i/'initiai) 
evolves in time as 

D 

|V^(t)) = ^c„e-‘^”*|i?n) (1) 

n—1 

with Cn = (n-l'i/initiai), where the Planck constant is set 
to unity throughout this paper. We assume that the 
energy gaps {En-Em,}i<Tn<n<D are all different 
This assumption holds true in nonintegrable hard-core 
Bose-Hubbard models (see below) or equivalent spin-1/2 
systems (see e.g., Ref. M)- We also assume that an 
effective dimension [l3| . which is the effective number 
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of the energy eigenstates involved in the initial state, is 
much greater than unity: 

Then the time-dependent expectation value of 
a few-body observable O, or ('(/;(<)|O|'0(<)) = 

{En\d\Ejn), is most of the time 
close to its infinite time average, since the off-diagonal 
(m ^ n) contributions cancel each other. The infinite 
time average is represented by the diagonal ensemble 
average 


D 

(0)de = ^ |c„|^On, (3) 

n—1 

where 0„ = is called an eigenstate expecta¬ 

tion value (EEV) and |c„p represents the energy distri¬ 
bution. Thus, Eq. ([3]) describes the physical quantities in 
the effective stationary state. We note that Eq. o ex¬ 
plicitly depends on the microscopic details of the initial 
state. 

The MCE gives a different average (0}mce for the ob¬ 
servable O, which is defined, with only two parameters, 
be., the central energy Ec and the energy width 6, by 

(O)mce = X! 

neI{Ec,S) 


where I{Ec, S) = {n\ En € [Ec — 6, Ec + ^]}, and the 
normalization factor Mec, 5 = 12n^i{Ec <5) ^ gives the 
number of the energy eigenstates in the energy window 
[Ec — S,Ec + ^]- In this paper, we consider only those 
cases for which (O)mce coincides with (0)de in the TDL. 
Then we may define the error of the MCE due to the 
finite-size effect by 

Error of MCE = (6)mce - (0}de (5) 


and numerically investigate this error in a concrete 
model. 

Model .— We use a one-dimensional model of HCBs 
with the nearest- and next-nearest-neighbor hopping and 
interactions. The Hamiltonian is given by 


L 

jjH -b blbi+i) + uhihi+i 

i=l 

L 

+ ^ [^—+ blbi+2) + , 

i=l 


( 6 ) 


where the periodic boundary conditions are imposed and 
hi is the annihilation (creation) operator of a HCB 
on site i with [bi,bj] = = 0 toi i ^ j, 

= 0 and {hi,hi} = 1, and hi = h\bi. The 
total number N of HCBs is conserved and chosen to be 
fV = L/3 in our numerical study. 


TABLE I. The number N of HCBs, the number L of the 
lattice sites, and the dimension D of the Hilbert space for the 
sector with the translational momentum P = 2tt /L which is 
used in our numerical study. 


N 

5 

6 

7 

8 

L 

15 

18 

21 

24 

D 

200 

1026 

5537 

30624 


For simplicity, the nearest- and next-nearest-neighbor 
hopping and the next-nearest-neighbor interaction en¬ 
ergies are set to unity [s^ in Eq. ([B]). In the follow¬ 
ing discussions, we consider quantum quenches by sud¬ 
denly changing the parameter u from 0 to it/. In these 
quenches, 77^“^ is nonintegrable and the energy level 
spacings obey the Wigner-Dyson statistics due to 
the next-nearest-neighbor contributions [ 3 ^ . 

Due to the translational invariance of our model, the 
Hilbert space is decomposed into L sectors labeled by the 
translational momentum P = 2'irmlL (m = 0,1,..., L — 
1), and we take the one with P = 27r/L, which can¬ 
not be decomposed into smaller sectors. The dimension 
D of this sector at each system size is shown in Ta¬ 
bleland the energy eigenstates in the sector are denoted 
by where the corresponding eigenenergies 

are arranged in an ascending order. All these 
eigenstates are obtained by numerically diagonalizing the 
Hamiltonian. 

Protocol of our numerical experiment. — We consider 
a quantum quench where the initial state IV'initiai) is an 
eigenstate \Eng) of The time evolution is governed 
by and an effective stationary state is eventually 

reached where the expectation value of a few-body ob¬ 
servable O is given by the diagonal ensemble average ® 
with c„ = {E^^'^\EnJ). Meanwhile, we calculate the 
MCE average of O ([1]) by setting S = 0.02L and Ec so 

that ('0i„itiai|ff^“^^|'0initiai) « (ff^“^^)MCE- Thus, we ob¬ 
tain the error of the MCE ([5|) for the given initial state 

We calculate the errors of the MCE starting from every 
eigenstate \Eng) whose “effective inverse temperature” 
/3„o falls in an interval [0,0.05]. Here and henceforth, 
the Boltzmann constant is set to unity. The effective 
inverse temperature /3„ of an eigenstate is defined 

by the equation En'^ = Z~^ J2m=i ’ where 

Z = . As shown in TableHIl the eigenstates 

thus chosen lie in the middle of the spectrum |38l | ; in this 
case chaotic states appear [ 3 ^ |4^ and thermalization 
occurs [H,Hi,mini. 

We investigate two local operators Oi = nih 2 and 
O 2 = fiifis representing the correlations in the numbers 
of HCBs between the nearest and next-nearest neighbors, 
respectively. We have confirmed that our results shown 
in the following are qualitatively unaltered for other two 
local operators, b\b 2 + H.c. and b\b 3 -|- H.c.. We note 
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TABLE II. The label n of the eigenstate that has the effective 
temperature closest to /3 = 0.05 and 0.00. The ratio n/D 
shows where each eigenstate lies in the spectrum. 


L _ 15 18 21 24 

/3 = 0.05 n 79 403 2140 11696 

_ n/D 0.40 0.39 0.39 0.38 

p = 0.00 n 96 485 2628 14578 

n/D 0.48 0.47 0.47 0.48 



{0)mce — (0)de 


(6)mCE — {d)DE 


FIG. 1. (Color Online) The distribution of the errors of the 
MCE (Eq. ((5|) obtained for each of the initial states in the 
quench of Uf = 0.4 for (a) [(c)] L = 18 (filled box) and 21 
(solid line) and (b) [(d)] L = 24 (solid line) for 6 = Oi [O 2 ]. 
The dashed curves in (b) and (d) are the least squares fits of 
the distributions with the Gaussian distributions. We note 
that the horizontal and vertical axes are scaled by the factors 
shown at the right-bottom and left-top corners of the panels, 
respectively. 


|(0)de — (0)mce| < 1 for O = Oi and O 2 because their 
operator norms are unity. 

The 1/D scaling of the accuracy .— Figured] illustrates 
the distribution of (0)de — (O)mce obtained for each ini¬ 
tial state in the quench of Uf = 0.4. FiguresHIa) and (b) 
show that, for Oi, the errors become markedly concen¬ 
trated near zero as the system size increases. Figuresdic) 
and (d) show that the distribution of the errors for O 2 
behaves similarly to that for Oi. We discuss, in the fol¬ 
lowing, how fast the width of the distribution vanishes 
by examining the root mean square (RMS) of the errors, 
which we call the accuracy of the MCE. 

The first main result of this paper is that, in the quench 
oluf = 0.4, the accuracy is proportional to 1/D for both 
Oi and O 2 as shown in Fig. |2]]a), where the accuracy is 
plotted against D with the error bars representing the 
estimation errors [i^. By conducting the least squares 
fits of the accuracy with a function f{D) = AD^, we 
find the exponent B to be — I.OOIq qj for O = Oi and 



FIG. 2. (Color Online) (a) Accuracy of the MCE for two local 
operators O = Oi (circles) and O 2 (triangles) at L = 15,18, 21 
and 24 (see Table [J for the corresponding D) with «/ = 0.4. 
The error bars show the estimation errors [1^. The data 
points are well fitted by the dotted line which has the slope 
of —1. (b) For Oi (circles) and O 2 (triangles, slightly shifted 
to the right for clarity), the exponent B obtained by the least 
squares fits of the accuracy with a function f{D) = AD^, 
where the error bars represent the 95% confidence intervals. 
The dash-dotted line indicates the value of —1/2 as predicted 
by the ETH. 


—O.OTIq;]]® for 0 = 02 with 95% confidence (—O.QQjlQ;^® 
for O = b\b 2 + H.c. and —for O = b\b 3 + H.c.), 
which clearly demonstrates the I/O scaling of the accu¬ 
racy. We have also conducted similar analyses for var¬ 
ious quench magnitudes and found the 1/D scaling for 
the range of 0.2 < Uf < 0.75 as shown in Fig. dKb). 
The discrepancies of the exponents from —1 seen for 
0 < Uf < 0.2 and Uf > 0.75 will be addressed later, 
and we here focus on the implications of the I/O scaling 
of the accuracy. 

We discuss the I/O scaling in terms of the number 
N of HCBs in our model with a general filling factor 
iz = N/L, which has been fixed to be 1/3 in the above 
discussions. The dimension O is approximately given by 
O Ri {^)/L Ri 10'>'('')^, where 7 ( 1 ^) = —v~^[i>\ogiQi> + 
(1 — ^) logiQ(l — v)] is a monotonically decreasing func¬ 
tion giving, in particular, 7(1/4) = 0.9768 • • • Ri 1. Thus 
the I/O scaling implies that, at 1/4 filling, the accuracy 
improves by one order of magnitude as we increase the 
number of HCBs by one. 

Two ETH upper bounds on the accuraey. — Before dis¬ 
cussing the underlying mechanism for the I/O scaling of 
the accuracy, we derive from the ETH two upper bounds 
on the accuracy of the MCE by assuming that the en¬ 
ergy distribution |c„p is localized, or its width is much 
smaller than the macroscopic energy scale such as the 
total energy. This assumption is needed for the state to 
be thermodynamically normal in the sense that the total 
energy is macroscopically definite. Then we point out 
that the upper bounds decrease proportionally to l/\/15, 
implying that the ETH alone cannot explain the 1/D 
scaling. 

First, we note that the EE Vs of a few-body observ¬ 

able O in nonintegrable systems are known to behave as 

On = / {En/L) -\- SO., 


(7) 
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where f{x) is a smooth function and 60 n represents ran¬ 
dom fluctuations around it [d^. Correspondingly, the 
error of the MCE is decomposed into two parts: 

(0)dE — (O)mCE = AOsys -I- AOrand, (8) 

where 

AOsys = |CriP/(E„/£) - /mCE, (9) 

n 

AOrand = ^ |Crap(50„ — JOmCE- (10) 

n 

with /mce = -^Ec^S^neiiEcS) fi^n/L) and 50 mce = 
•^Ec,S SnG/(-Ec.<5) 

Second, we note that AOgys is negligible if the energy 
distribution |c„p is sufficiently localized. In fact, under 
this condition, the Taylor expansion of f{En/L) up to the 
first order is sufficient and we have AOgys ~ 0 because we 

choose Ec so that (■i/'initial|l?^“''^ l^/’initial) ~ (ff^“^^)MCE- 
We discuss the higher order contributions in Discussion 
below. 

Thus, if the energy distribution |c„p is sufficiently 
localized, the error of the MCE is dominated by 
AOrandi which is bounded from above strictly by Aq = 
2max„ |i50„| and roughly by the standard deviation (jq 
of {50n}n- We note that Aq and ao should be calcu¬ 
lated in an energy window where Icnp’s are significantly 
weighted instead of the entire spectrum. In the following 
discussion, we take the window to be the microcanonical 
one because the width |c„p has turned out to be smaller 
than 5 ioY Uf < 0.75, where we discuss the 1/D scal¬ 
ing. We call both Ao and cro the ETH upper bounds 
because they are commonly used as the indicators of the 
ETH in the strong and weak senses, respectively, which 
imply Ao —>■ 0 and cro —0 in the TDL |27l - [^ H . 1^ . 

It has recently been shown that cro approaches 
zero proportionally to I/'/D in nonintegrable spin sys¬ 
tems [13, m, . We have also obtained the 1 / 'Jd scal¬ 
ing for both ETH indicators Ao and uo as illustrated in 
Fig. [3] These results imply that the ETH alone cannot 
explain the l/D scaling found in our numerical study. 

The underlying mechanism for the l/D scaling. — The 
second main result is that yet another indicator 


(To = 


a/ deff 


( 11 ) 


can describe the accuracy as illustrated in Fig. [31 This 
indicator involves, in addition to ao, an extra suppres¬ 
sion factor l/y/des- Since des represents the effective 
number of nonzero terms on the RHS of Eq. (HOD, the ex¬ 
tra suppression factor 1 / y/dfs implies that there is little 
correlation between the terms. 

It has been proposed as an alternative to the ETH [13, 
[H, [13, [i3| that the absence of correlation between the 
energy distribution |c„p and the EEV 0„ suppresses the 
error of the MCE. Indeed this is the case for our setup 
since the initial states are chosen independently of the 




FIG. 3. (Color Online) Accuracy of the MCE (circles), strong 
(filled squares, Ao) and weak (open squares, ao) ETH indica¬ 
tors, the indicator for the no-correlation model (triangles, do 
defined in Eq. (1111) 1 for (a) Oi and (b) 62 with Uf — 0.4. The 
solid and dashed lines with slopes —1 and — 1 / 2 , respectively, 
are guides to the eye. The numerically obtained accuracy is 
best fitted by the no-correlation model. 


observables Oi and O 2 . However, this mechanism alone 
cannot explain why the accuracy is so good as Eq. m 
if |c„p’s are correlated with each other. 

Let us derive Eq. m by introducing the “no¬ 
correlation model” in which there are no correlations 
between |c„p’s except for the constraints J2n = 1 
and ~ 1/c^eff- To be more specific, it is given 

by |cnP = gnXn, where qn represents a smooth pro¬ 
file of |c„|2 satisfying En 

and = 0 for n ^ I{Ec,S), and {xn}n=i posi¬ 
tive random variables whose mean and standard devi¬ 
ation are both unity. These conditions are known to 
hold in the middle of the s pec trum in chaotic systems 
with Qn being Gaussian [ill l48l. 1^ . The no-correlation 
model leads to (AOrand)^ = ((To)^/(ieff+(^0‘^-^OMCE)^, 
where SO‘> = dnSOn, and (cr^)^ = J2n ^n(^On)^ with 
fn = Here ttt denotes the statistical aver¬ 

age with the probability distribution P{xi,X 2 ,... ,xd) = 
p{xi)p{x 2 ) ■. .p{xd)- Since {g„}n and {rn}n are smooth 
and normalized profiles within the microcanonical win¬ 
dow, we assume <50^ = (^Omce and (ctq)^ = (Tq. Then, 
we obtain (AOrand)^ = (To/(ieff, which implies that 
Eq. (fTTl) gives an estimation for A Grand ■ 

Our result that the no-correlation model explains the 
l/D scaling of the accuracy implies that we cannot 
induce correlations between jcnp’s through a sudden 
change of a single parameter in the middle of the spec¬ 
trum where many-body energy eigenstates show chaotic 
behavior. While we can control the total energy by 
changing the parameter, we cannot manipulate indi¬ 
vidual many-body eigenstates whose landscape changes 
drastically from one to the neighboring one for noninte¬ 
grable systems. 

Discussions .— First, we discuss how the l/D scaling 
is modified in larger system sizes. As we increase L fur¬ 
ther, there might exist an Lupper at which AOgys be¬ 
comes comparable with AOrand and can no longer be 
ignored. This is because, whereas AOrand decays as 
l/D, AOgys decays only algebraically with L for ther- 
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FIG. 4. (Color Online) Accuracy of the MCE (circles) and 
indicators for the no-correlation model (triangles, ao defined 
in Eq. HU) and the weak ETH (open squares, uo) for Oi 
with (a) u/ = 1 and (b) uj = 0.05. The solid and dashed 
lines with slopes —1 and —1/2, respectively, are guides to the 
eye. 


modynamically normal states as follows. By making 
the Taylor expansion of f{En/L), we have AOgys 
2~^ f" {Ec/L)[AE^^ — AEIiq^]/L'^ because the zeroth- 
and first-order contributions cancel out due to the nor¬ 
malization condition and the fact that Eq is chosen to be 
equal to ('(/'initial IV'initiai). Here AEl,^ and AEI^^^ 

denote those the energy fluctuations in the DE and MCE, 
respectively, that are proportional to L, and we obtain 
AOsys oc L~^. Even if we perform the fine-tuning of 5 so 
that AE^g = AE^Qg, the contributions from the Taylor 
expansion of f{En/L) at all orders cannot be canceled in 
general, and AOsys decreases only algebraically with L. 

The Tupper should be smaller when the energy distri¬ 
bution |c„p is less localized and the higher-order con¬ 
tributions in the Taylor expansion of f{En/L) become 
important. In fact, for the large quench with m/ = 1, the 
deviation from the \/D scaling is seen at L = 21 and 24 
as shown in Fig. ITJa). The deviation may be interpreted 
to be the cause for the deviation from the 1 /D scaling at 
the largest quench in Fig. [5Db). 

Second, we discuss how the 1/D scaling changes in 
smaller system sizes. As we decrease the system size 
with the quench magnitude held fixed, the quench en¬ 
ergy becomes smaller than the energy level spacings and 
we enter the regime where deg 1. Thus, for a given 
quench magnitude, there exists Tiower below which the 
1/D scaling disappears. In this case, Eq. (ITTI) reduces 
to the indicator of the ETH in the weak sense and the 
accuracy is proportional to l/y/D. This crossover be¬ 


tween the 1/D scaling and the I/'/D scaling is seen for 
M/ = 0.05 as shown in Fig. SDb) [H^. Thus, the lack 
of mixing of numerous energy eigenstates is the cause 
for the crossover from the 1/D to I/'/D scalings for the 
small quenches as can be seen for u/ < 0.2 in Fig. [IDb). 

Thus, we find three regimes of the system size for a 
given quench magnitude: (i) L < Tiower where the accu¬ 
racy is described by the ETH and proportional to I/'/D, 
(ii) Flower < L < Fupper where the accuracy is propor¬ 
tional to 1/ZI due to the absence of correlations between 
the many-body eigenstates, and (iii) Fupper < L where 
the accuracy improves only algebraically with L because 
AOsys rather than AOrand dominates in the error of the 
MCE (see Eq. (H])). 

Our findings imply that there exists an exponential en¬ 
hancement of the accuracy in the small system sizes and 
the MCE can describe equilibrium states very accurately 
even in small isolated quantum systems. In fact, the ac¬ 
curacy of the MCE reaches or even better in the 

system with only 8 HCBs on 24 sites as shown in Figs. [3] 
andlH 

Conclusions .— We have numerically investigated the 
accuracy of the MCE in interaction quenches for a nonin- 
tegrable hard-core Bose-Hubbard model. We have found 
a regime where the accuracy improves proportionally to 
1/D (see Fig. [2|). This rapid improvement of the accuracy 
implies that quenching a single parameter cannot induce 
correlations between the numerous many-body eigen¬ 
states since they depend nontrivially on the parameter in 
nonintegrable systems. As we increase the system size, 
there are three regimes where the accuracy scales as (i) 
I/'/d (ETH regime), ( ii) 1/D (no-correlation regime), 
and (iii) algebraically with L. Due to the regimes (i) and 
(ii), where the accuracy improves exponentially with L, 
the MCE can describe the equilibrium states quite accu¬ 
rately even in small systems. 
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